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ABSTRACT 


The work of Lewis and Walecka on the structure of the giant reso- 
nance in Carbon 12 and Oxygen 16 has seemed to indicate that the 
single-particle, shell-model theory gives better results than collective 
models -- in particular the Goldhaber - Teller and Steinwedel - Jensen 
models. The transverse form factor as a function of momentum 
transfer is observed experimen elle’ fall to a minimum and rise 
again in the range from zero to 0. 6 Se » his dip 1s predicted by tic 
Brown particle-hole theory model but 1s not predicted by either of these 
collective models;* The purpose of this paper is to develop a collective 
model which could possibly give this dip in the form factor, Further 
work is necessary to determine if this model can actually do so. A 
semiclassical treatment is used to find the matrix elements for the 
assumed model. The transverse form factor is then calculated for two 
nuclear charge density functions -- constant density and the Fermi 


distribution, 
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1, Introduction. 

Electron scattering is a very useful tool for the study of nuclear 
structure, since the cross-section for electron scattering from a 
nucleus can be related to the basic geometric details of the structure 
of the nucleus. The incoming electron is assumed to interact with the 
nucleus by means of the electromagnetic field. This interaction is 
known, whereas the nuclear interaction (strong force) is not. 

the parameters of electron scattering are the imitial momentum k 
and initial energy E of the incoming electron, which are controlled by 
the experimenter; and the final momentum k‘’ and the final energy Be 
which are measured. The cross-section can be found as a function of 
the momentum transfer q = k - k’ from the electron to the nucleus 
and of the excitation energy, which can be found from E, BE’, the 
scattering angle, and considerations of the nuclear recoil. One can then 
Vaty q for a fixed excitation Energy (i. ¢,, for an excitation of the 
nucleus to a particular level), and get the cross-section as a function 
of G (where q's | q | ) for that level. This cross-section is related to 
the Fourier transform of the charge and current densities of the nucleus, 
Thus, we can get information about the charge and current densities of 
the nucheus from electron scattering. 

A so-called ''dipole state'' is a feature of the electron scattering 


cross-section for most nuclei. The dipole state is a state strongly 


excited by photons of the appropriate energy (usually 10-25 Mev). It 
shows up as the ''giant dipole resonance"', which is a resonance peak 
in the cross-section as a function of energy, The dipole state is so 
named because it has a large electric dipole matrix element with the 
ground state, As a result of this last fact, these states lead to strong 
photon absorption cross-sections which are orders of magnitude larger 
than those of other multipoles, 

lf nuclear recoil is 16en0recd; oe cross-section for electron 


scattering can be written in the general form [2]: 





2 a 
dg 8a k’ + ; 2 2 
ii ; 7 | ji 2kk’ cos (6/2) fo (q) 
q q 
Zick A 2 2 2 2 2 
5- sin (0/2) [q + 2kk" cos (@/ 2) J f.  (q) | 


where q is the four-momentum transfer and @ is the angle between the 
incoming and final electron momenta. fe and f are the longitudinal 
and transverse form factors, ~Recardless of the model used, the cross - 
section will have this form and will depend on these two form factors, 

If nuclear recoil is taken into account, the cross-section is simply 
multiplied by a factor 


] 
k’ -k cos Q 


ee ow 
The transverse form factor can be easily separated out, since for 


O 
8 = 180 , the coefficient of the longitudinal form factor goes to zero. 


Figure | shows the experimental results for the square of the 
; EZ 

transverse form factor for electron scattering from CC as a function 
of momentum transfer, along with the predictions of two collective 
models (the Goldhaber - Teller and Steinwedel - Jensen models) and two 
particle-hole shell models (the Brown and Gillet - Melkanoff models), 
Figure 2 shows the experimental results for the transverse form factor 

16 cases 
for 0 and the predictions of the Brown and Goldhaber - Teller models, 


For both e= and glo 


, the observed behavior is more closely predicted 
by the shell models than by either of these two collective models, Lewis 
and Walecka assert that this is conclusive proof that the collective mode 
description of these excited states is not adequate, and that the particle- 
hole shell model description is correct. The purpose of this paper is 
to develop a collective model which can predict the observed behavior 
of the form factor. 

In the Goldhaber - Teller model, the neutrons of the excited nucleus 
arc assumed to move against the protons as a unit, thus giving the 
effect of an harmonic oscillator, with the restoring force being due to 


changes in the overlap of the protons and the neutrons [cj]. The 


interaction Hamiltonian is assumed to be 


— 
I 
_— 


= 5 Pi] Q° 
= 5 to , 


where w is the frequency of oscillation, u= NZM/A is the reduced 
mass of the system, Q is the coordinate representing the separation 


of the center of mass of the proton system from the center of mass of 


the neutron system, and P is the canonical momentum associated with 
the coordinate Q. PandQare quantized and the displacements are 
assumed to be small. The charge and current densities are found from 
the motion of the protons. 

In the shell model theories, a single nucleon is assumed to be 
excited {roma closed sicll by the interaction with the incoming clect ron: 
This mucléonvthenimteracts with the -resultmp- “hole” leit in tie sheulat 
previously occupred,. “rhe nucleon and tne hole are treated as point 
particles. The excited state of the nucleus is assumed to be a linear 
cOmipindlion of the single-particle excitations, and the ¢nar pe, current, 
and magnetization matrix elements between the single-particle states 
are found, 

In this paper, a collective mode is assumed and a semiclassical 
treatment is used to find the matrix elements of the Hamiltonian de-~- 
scribing the interaction between tie clectron and the nucleus. ihe 
cross-section is then found to depend on the two expected form factors, 
These form factors are then calculated for a nucleus of constant density 


and for a nucleus whose charge distribution is the Fermi distribution. 
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Zz. TiemModel,. 

The collective mode assumed is an harmonic oscillation of the 
protons and neutrons in the @-direction about some axis, so that the 
Charge demsity of the excited part of the nucleus is: 


o(r) = Ze f(r) (1 + a cos@ eee ) 


where Z’ is the number of protons involved in the excitation, and 

Z‘e f(r) is the ground state charge density of this portion of the nucleus, 
It is assumed that the neutrons! motion is equal and opposite to that of 
the protons, so that the net angular momentum is zero, At this point, 

no specification 1s rade as to whether the entire nucleus or a portion 
only is excited, a]though it may later be shown that in order to get the 
desired results for the transverse form factor, it is necessary to assume 
that only an outer core: is excited, with an inner core of nucleons re- 
maining inert. w, the frequency of the oscillation, is equal to the 
encroy transicrred to the:nucléeus by theanteracting e€lecitron. ‘Ihe 
system of units employed inthis paper is that in which c=T =1. 


From charge conservation, 


—= — op 
V : J _ = e e 
The assumed mode has 
1 =) are 
96 


where es is the unit vector in the @ direction, Thus, the charge 
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conservation equation is: 


Lw t 
—— — (J Sin6)he= -iwaZ’e f(r) cos® Pe 


This leads to: 
Lat 


et) = -5 igage = f(r) sin@ e 


Thercaleulation of the matrix elements A. is done semi-classically, 


The matrix element is (5 | 
3 
H. = | 3,094, 69 dx . 
fi U UL 


It is assumed that the nucleus is scattered by the electric field of the 


electron, which can be written as 





where q is the four-momentum transfer and uy and u, are the final 
and initial plane-wave electron states[5]. Vis the normalization 
volume for these states, It is further assumed that J (x) for the 

U 


nucleus can be written as just the classical charge and current densities 


for the model, Thus: 





S 4fJ1i oles - 3 
Hie. = [J (xa xo xo ae gee J x) u. Y wuld. 
a b q V be Big 
Lwt 
But J (x) is just j (x) e” , and 
4 + 
| io oe aod x aoe . 


i2Z 


so what must be calculated is the Fourier transform of the current 


density: 
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3. Caleulation of Fourier Transforms. 

There will be three normal modes for this model: one for which 
the axis of vibration is parallel to q (the longitudinal mode) and two for 
which the axes of vibration are perpendicular to q (the transverse modes), 


For the longitudinal mode, let the z-axis denote the axis of vibration. 


Then: 
S Zz 
j_(a) = ae: iN Stn Gt deociemeO some searches 
=O 
Likewise: 
2 
jy (a) = c} an f(r) sin §@ cos®@ sing EID acl ae 


= 0 


Thus, for the longitudinal mode, both transverse terms are zero. 





However: 
; 5 Oye. l 
j, (a) = - c | eS itr) ssa 9 Fe Cll 
8 
—— nS | £() PSsinqrcy=-qr il cos (qr jiijwdr 
q 
- ATT; / 
TliaZ ew f, (q) 
where 
] 
fF (a) ears | Csinax - qr cos(qr)]dr. 
q 


The fourth component of j (q) is: 
U 


9 ; 
ia | £(r) r sin@ cos@ Sp ha na 


/ 
-4MaZ’eq £, (q) 


i [p(w oes oe 


It 
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Thus, for the lonyitudinal mode, 
j (q) = 4MiaZ’ew f, (q) (0,0,1, 4). 
U L Ww 


For the trunsverse modes, the calculations are done in the coordi- 


nate system in which the z-axis is the axis of vibration and the y-axis 





1s Iptine direction of q. For one of these modes: 
| 2 
1 p (q) La fir) sin@ cos®6@ et” a drd@d¢ = O 
- I Le | iqr sing sing 
ean) —a © | f(r) Sin @ cos@ cos¢ e drd@d¢ 
0 
2 ; 
10a) : c | frye sin 8 cos@ sing eI eee drd@d¢ 
0 
5 yeas iqr sin@ sing 
1, (4) = mom ican sin @ ec drd@d¢ 
J 
- ne | f(r) [qr - sin(qr)] dr = - 2[liaw eta) 
q 
where 
l : | 
fq) | | fe) [ sin (qr) -qr]dr. 
q 


For the modes in which the axis of vibration is perpendicular to q , 
the only nonzero Fourier transforms are those in the directions of 
the axes of vibration. The overall current density Fourier transform 


can be written: 
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where F 1s the polarization axis for thesexcited mucous ps2. jon >, 


with the 3-axis being the axis parallel to q mee 
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4. 


Calculation of Cross-Section, 


The interaction Hamiltonian is now: 





Z 
161] 5 yn Pp - 
= = 2 ft (q) ee N u. 
Lb 
q V 


The cross-section depends on the square of the Hamiltonian, 


averaged over initial electron spin states and summed over final 


electron spins and nuclear polarizations, 





This 1s: 
Z l fem Zoe" w P - é 
jH|" = 7 f Ge rene \f (q)u,y u. | 
i, f, P q Vv Ih b 


The summations lead to a trace calculation, as follows; 
Pp = os Pp se 
{ i = 
ae Ga "6 a, te Way » 


: ‘ P* P | 
rr | Cay, y, Vga up ) fy f : 


UL 


ei toweexauiine the Quantity in parentheses, Ifthe vs and 


Vg were Ia versed, this would just be 


Tr Uy Gy uy), 


which is [5]: 


Ly 


] / / a, 2 
mam—— ([kk’ + kk’ - (k *k’ +m )6 jj. 
EE? wb Vv Vou = = Lb V 


From the commutation relations for the gamma matrices; 


et Yy Y4 ue) 





ne : 
Se en epee eo gerry 
EE wu Vv Vou _- — Uv ie ee 


The form factor matrix is: 


0 0 9) 0 
PS0*:0- 0 0 0 O 0 
Lee Z 
i, ras (q) Oo Oo a (q) 0 0 l iq/w 
00 0 0 
228 
0000 00 -ig/w gq /w 
aes t 
=- —f A i A 
zt (a) vin i !w on 


There are two traces which must be evaluated, These are: 


The transverse coefficient: 





t 2 2kk ‘ 2 
Tr(B A a 4 sin (0/2) E - a cos (0/2) | 
Vow 2 
q 
The longitudinal coefficient: 
Z 
iu oa 2 
Tr{B A = ys cos (9/2) 


V Vv 
bo Wb 2 


where @ is the angle between the incoming electron momentum k and 
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the outgoing electron momentum k‘, The square of the Hamiltonian, 


averaged and summed over states, iS: 


~~ a 
2 16M aZ‘e 2 eee ne 
ee rere) w Cree { 2kk’ sin’ (@/2) 
q V 





x | | ‘ ans oe (9/2) | i (q) 


fs 
q 
4 
+ 4kk" q 2 2 
<P (6/2) 7 (a) 
q w 


The amplitude a of the oscillation can be calculated (classically) 
from eneryy considerations, At the time during the oscillation such 
that the charge density of the nucleus is its ground-state configuration 
(wt (n + >) ), aii the enerpy transterred to the nucleus is inthe 
kinetic cnergy of the neutrons and protons, The velocity distribution 
ol ithe Cxcited protons is 


vir) = j(r)/Z’ef(r). 


war sing. The mass density of the excited 


NM] 


Tit@iratiov1s vir) = 
protons is Z’M f(r). Thus, the kinetic energy distribution of the 


excited protons 1S: 


Z an on 2 
= MZ‘ é(r) ae on ae “ Z'M f(r) w ar Sin @ 
Since half of the transferred energy will be in the protons (with the 


other half being in the kinetic energy of the neutrons): 
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pee 
2 2 Z 
5 _ MZ 0 ae = a pir). + sin @ r sinO@Q d@Q@dr d¢ 
This gives 
zZ 3 
a —_—_—_—_—oo ’ 
2T1 Z‘Mul 
4 
where I = jf) Pv dit.. 


The cross-section is, from the well-known Golden Rule of Ferm1, 


2 1 

sess, Se: ee gee) Gee 
final 
states 


The flux is just k/ VE. . The summation over final states is just 


3 : Doe 
an integral over d k" for the case at hand, The 4-function is just 


Sai eae =e) ae hus: 








Wea, ee 
V 2 2 
do = ponies —— |H,|° 6(k-k’-0) k’% dk’ do 
(211) 
2 
VE, k’ 5 
= ——>— |[H,.|° do , 
Z| aoe 
Oh iame fi k - k Ww 
, | Ls , 
Finally, putting psa into this equation, we obtain: 
d gone! Z 
oT ar, SS | 2kk’ any (9/2) ( [ee cos (9/2) | 
qk 20 Z2'MI 
Bee 
q 2 2 
+ nea aes to COS (9/2) i 
q W 
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Note that these form factors must be multiplied by a constant 


term 


W 


Fa 
21° Z'MI 
in order to compare them with the form factors shown in Figures 1] 


and 2. 


el 


5. Uniform Snell Wviodel: 
If the nuclear charge density is assumed to be uniform, with the 
nucleons within a shell of inner radius ry and outer radius r5 being 


excited by the interaction, we have: 


f(a) = 3 | sin (ar) - ar, cos (qr. ) ~ [ sin{qr, ) - qr cos (ar) ] | 
l 
i (qj). = a [ sin(ar,) ~ qr, - [sin{qr, } ~ ar] | : 


q 


Fireure-3-shows the resultsior the squareof the transverse torm 
16 . 
factor versus momentum transfer for Q , with the assumption that the 
entire nucleus is excited (r) = 0), and the same form factor with the 
: , ee 4 
assumption that an inner core consisting of a He nucleus remains 
inert. ,Dhewtipures wsed ior the Outer radii for these calculations were 
16 ; 4 

3.41 fermi for 0 and 2.09 termi forHe | 3 |. 

If one aSsumes an even thinner shell is excited, the form factor is 
reduced by orders of magnitude and falls off to zero with increasing q 
even faster, Figure 4 shows the results for a shell of thickness 0. 1 
fermi. No matter how thin the shell, the uniform shell model cannot 


give the observed form factor behavior with q for a reasonable outer 


radius. 


Ze 


6. The Fermi Distribution Model. 


Now Ict us assume that the charge density function is just the 


Mermn) distribution: 


Cc 


a 


leet? . pete 


The integral can be done by breaking it up into two parts and writing 


ach part as an infinite series, as follows: 





im 2 ' 
f{ (q) Pirate) - qrdr = | t lhe | f(r) 
{ 
O 0 
X [sin(qr) - qr ] dr 
EOP eo 
SEN : P as na (r -a) 
| + € O 
Otel 
l = nt] -na(r-a 
- (r on _ Le ( - 1) € ( 
y oe 
Ui eresultior the term factor rs: 
2 2 2 sin (a) 
a ] a hi 2 
(q) = ae 5} a) - Z 
q ba q 
aes ? 2 cos qa 
oo ac Z a 
n 7 
a | Soe | 
| q +8 @ 


For excitation of a shell, the form factor is 


where a is the outer radius of the shell and b is the inner radius. 
The square of the transverse form factor for oe is shown in Figure 5, 
using a = 3,64 and these combinations for banda: b = 3.0, 
ge-=e 25. ewoeo, woe 1 end! =.5--0 5.0 —2c0le acai theta mem 
the shell, the smaller the form factor and the faster it drops off with 
increasing gq. “For large-a.. the square of the form iactor 16 smaller 
and does not drop offas fast. In fact, for very large aw, one gets just 
the Uniform shell. “Note that the * curve for a = -cO ‘actually i¢uosses 
the w@ = | curve, This leads one to suspect that if one asstimes an outer 
boundary for the excited shell with, say, qa = 1 and an inner boundary 
with very large a, for a thin enough shell one could get a curve for 
pe | G which has the desired minimum at about “q =,0.435. 

Figure 6 shows the results for three trials of this omposite model. 
Hor a thick shell. (a =-3,.64 and) b =-2, 0) swithrd = 1.° the form factor is 
larger thanctor the same shellywith sy —-110,- “lhevform ta ctor becomes 
larger for a thinner shell. However, for the combinations of a, b, anda 


used here, the form factor doés not get the desired behavior, 
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7. Conclusions, 

None of the assumed charge distributions gives the desired behavior 
for the form factor. This model gives somewhat the same results as 
the Goldhaber - Teller model; that is, the function | f(q) | C is a 
motonically decreasing function of q within the range of interest. 
Further work is necessary on the composite model before this model 
can be completely discarded, however. For instance, if one assumes 
the nucleus as a whole has a Fermi charge distribution and that the inert 
core is 4 uniform shell with a radius equal to the radius of the half- 
height of the Fermi distribution, what one has is essentially a classical 
Precatimentvol the particle-hole theory, Figure 7 1)lustrates tne 
"particle'' and the "hole'' associated with this model. This idea should 


be pursued further. There is a strong possibility that it will give re- 


sults similar to those of the successful Brown particle-hole theory, 
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